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Abstract 

A detailed study is made of the effects of variations in 
lamination and material parameters of thin-walled composite 
frames on their vibrational characteristics. The structures 
considered are semicircular thin- walled frames with I and J 
sections. The flanges and webs of the frames are modeled by 
using two-dimensional shell and plate finite elements. A mixed 
formulation is used with the fundamental unknowns consisting 
of both the generalized displacements and stress resultants in 
the frame. The frequencies and nxxles predicted by the two- 
dimensional finite element model are compared with those 
obtained from experiments, as well as with the predictions of a 
one-dimensional thin-walled beam finite element model. A 
detailed study is made of the sensitivity of the vibrational 
response to variations in the fiber orientation, material proper- 
ties of the individual layers, and boundary conditions. 

Introduction 

The physical understanding and the numerical simulation 
of the dynamic response of laminated anisotropic structures 
have recently become the focus of intense efforts because of 
the expanded use of fibrous composites in aerospace, automo- 
tive, shipbuilding, and other industries, and the need to estab- 
lish the practical limits of the dynamic load-carrying capability 
of structures made from these materials. References 1-45 are 
indicative of the general interest and efforts focused on various 
aspects of the vibration of structures. Experimental studies 
have Ixren performed on the free vibration and impact response 
of thin- walled composite frames and stiffeners. 1 4 One- 
dimensional theories have been developed for the static, 
vibration and buckling analyses of thin-walled frame 
structures/ K However, no systematic assessment has been 
made of the range of validity of the basic assumptions of these 
theories. Approximate analytical and numerical techniques 
have been applied to the study of the vibrational response of 
isotropic and composite stiffeners.'* 21 Only a few 
publications 22 24 examine the effects of variations in lamination 
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and geometric parameters of composite panels on their vibra- 
tional characteristics, and these publications do not consider 
thin-walled composite frames. 

The present study summarizes the results of a recent 
study (Noor, Carden, Peters) 24 on the effects of variations in 
the lamination and geometric parameters of thin walled com- 
posite frames on their vibrational characteristics (frequencies, 
and energy components associated with different modes). The 
frames considered are semicircular, made of thin- walled 
graphite-epoxy material with I and J sections and have a 
36-inch radius (see Fig. I ). 

Analysis 

Computational Models 

Two computational mtxlels are used for the thin-walled 
composite frames study. In the first model, the flanges and 
web are modeled using two-dimensional shell and plate finite 
elements. The second irnxlel is a finite element discretization 
of the one-dimensional Vlasov's type thin-walled beam theory. 
Henceforth, the two models will be referred to as two- 
dimensional (2D) and one-dimensional (ID) finite element 
models, respectively. 

Mathematical Formulation 

a) Two-dimensional models . The analytical formulation 
for the two-dimensional nxxlels is based on the Sanders- 
Budiansky shell theory with the effects of transverse shear 
deformation, and laminated anisotropic material response 
included. A mixed formulation is used with the fundamental 
unknowns consisting of both the generalized displacements and 
the stress resultants in the frame (see Fig. 2 for the sign 
convention). 

Bicubic shape functions are used to approximate each of 
the generalized displacements and the stress resultants, The 
number of displacement nodes in each element is 16. The 
stress resultants are allowed to be discontinuous at interelement 
boundaries. The total number of stress-resultant parameters in 
each element is 128. The element characteristic arrays are 
obtained by using the two-field Hellinger-Reissner mixed 
variational principle. 
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b> One-dimensional models . The analytical formulation 
for one-dimensional models is based on a form of Vlasov's 
thin- walled beam theory with the effects of flexural-torsional 
coupling:, transverse shear deformation, and rotary inertia 
included. The fundamental unknowns consist of seven internal 
forces and seven generalized displacements of the beam (see 
l ig. } for the sign convention). The element characteristic 
arrays are obtained by using a modified form of the llellinger- 
Reissnei mixed variational principle. The modification con- 
sists of augmenting the functional of that principle by two 
terms: I ) the Lagrange multiplier associated with the con- 

straint condition relating the rotation of the cross section and 
the twist degrees of freedom; and 2) a regularization term that 
is quadratic in the Lagrange multiplier. Only C' continuity is 
required for the generalized displacements. Lagrnngian inter- 
polation functions are used for approximating each of the 
generalized displacements, internal forces and Lagrange 
multiplier. The polynomial functions for the internal forces 
and the Lagrange multiplier are one degree lower than those of 
the generalized displacements. In the present study quadratic 
polynomials are used in approximating the generalized dis- 
placements. Linear polynomials are used in approximating 
each of the internal forces and the Lagrange multiplier. The 
internal forces and the Lagrange multiplier are allowed to be 
discontinuous at interelement boundaries. Lor each element the 
total number of generalized displacement parameters is 21, the 
total number of internal force parameters is 14 and the total 
ntimlvr of Lagrange multiplier parameters is 2. Noor, Peters 
ami Min^ present the fundamental equations of the thin- walled 
beam theory. 
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Lor quasi-isotropic laminated composites, numerical* 
experiments to be described subsequently have demonstrated 
that reasonably accurate results can lie obtained using the 
one -dimensional mode! when the laminated composite is 
replaced by an equivalent isotropic material with the following 
Voting's and shear moduli: 


Anderson 2 *, and Noor and Peters 27 for the two-dimensional 
models, and Noor, et al 2 * for the one-dimensional model. 

Vibrational Sensitivity to Variations in Lamination and Mate- 
rial Parameters 


1’he expressions for the sensitivity derivatives of the 
frequency and response vector with respect to the lamination 
and material parameters, Xj, of the composite frames are given 
by: 2 * 


elements 
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and cv are multipliers given by: 
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In Fa|s. 4 to 7, the eigenvectors are assumed to be nor 
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malized with respect to |M|. i.c. 
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The expressions for the total complementary strain 
energy of the frame. ll r , and its derivatives with respect to Xj, 
are given by: 


where A )( and A n are the extensinnal stiffness in the x direc- 
tion. and the in-plane shear stiffness used in the classical 
lamination theory, respectively: and h is the total wall thickness 
(ol the flange or web). This approximation was adopted in the 
present study. 

Finite I dement Equations 
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1 he finite element equations for each individual element 
ot the II) and 2D models can lie cast in the following compact 
form: 


(|K| -w J |M|){Z| =0 (.1) 

where |Z) is the vector of the element degrees of freedom; co is 

♦ * 

the frequency of vibration; |K| and |M| are the generalized 
stillness and mass matrices. The explicit forms of matrix 

♦ t 

arrays associated wiili |Z|, |K| and |M| are given in Noor ami 


For llie purpose of obtaining analytic derivatives with respect 
to some of the lamination parameters, such as the filler orienta- 
tion angle of different layers, it is convenient to express ^ in 

. , dlFI 1 r „ 

terms of as follows: 
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The matrix is evaluated using the analytical derivatives 

of the material stiffness matrix of each laminate (flanges and 
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web). ’Hie material stiffness matrix of the laminate is given in 
Jones. 20 

Experimental and Numerical Studies 
Apparatus and Test Procedure 

a. Sped mens . Two specimens, shown in Fig. 4, were 
tested, an I -sect ion and a J-section frame. Nominal dimensions 
of each cross section are shown in Fig. I. Weight of the frame 
sections was 3. 181 and 4.085 lb. ( 1.443 and 1.853 kg) for the 1 
and J frames, respectively. The frame sections were made from 
AS4/5208 graphite/epoxy unidirectional tape layed up in a 
manner which resulted in essentially uniform stiffness proper- 
ties in the circumferential direction (i.e., the stiffness coeffi- 
cients are independent of 0). The material properties for the 
individual layers are given in Fig. I. The laminate stacking 
sequence for the l-section was |±45/0/90|, and |±45A)/9()| 2 , for 
the J-section. Each frame section was semicircular with a 
diameter of 72 inches (1.8288 m.). Bonded to the outside 
flange of each frame was a sixteen- ply | ±45/0/90 | 2s , quasi - 
isotropic skin made of the same material. The frame sections 
were constructed so that the skin would extend 0.5 inches 
(0.0127 m.) beyond each side of the bottom flange of the 
frame. Measured dimensions were used in one of the finite- 
element models and results were compared to nominal dimen- 
sion results and the experimental data. 

b. instrumentation and test method . A photograph of 
the lest equipment and composite frame specimens is shown in 
Fig. 4. The ends of the frame sections were potted in a fixture 
which was boiled ton large steel beam backstop. 

An air-shaker, connected to an air compressor, was used 
to excite all test specimens. Excitation was both in-plane 
(radially), and out-of-plane. For in-plane excitation, the shaker 
was positioned so that the pulses of air struck approximately 
along a normal to the surface of the skin. For out-of-plane 
excitation, a piece of styrofoam was attached to the side of the 
frame by double-sided adhesive tape. Pulses of air struck the 
flat face of the styrofoam along a normal to the face. The 
position of the air-shaker was adjusted if the excitation was 
striking on a node. 

A miniature accelerometer was attached at a fixed 
location to the frame sections with double-sided adhesive lape. 
Output from the accelerometer was amplified and displayed 
along the vertical axis of an oscilloscope. Natural modes were 
determined by tuning the excitation frequency of the air-shaker 
to produce an acceleration maximum on the vertical deflection 
on the oscilloscope. Output also passed through a low pass 
filter and was displayed as vibrational frequency on a fre- 
quency counter. ' 

A handheld velocity probe was moved along the frame to 
determine mxle locations and mode shapes. The output of the 
probe was displayed along the horizontal axis of the oscillo- 
scope. 'Hie probe and accelerometer outputs combined to 
create a Lissajous pattern on the oscilloscope, A phase shift in 
the Lissajous pattern occurred when the velocity probe passed 
over a node. 


Since manual equipment was used in mapping the nodal 
locations during the vibration survey of the frames, only nodal 
lines associated with gross in-plane, and gross out-of-plnne 
motions were monitored. Other nodal lines, associated with 
localized deformation patterns were not surveyed. These 
localized deformations were noticeable in some of the higher 
vibration modes, with complex deformation patterns and/or 
strong coupling between in-plane and out-of-plnne motions. 

Finite-Element Grids 

Two-dimensional models were generated for the actual 
frames (test specimens) described in the preceding subsection, 
as well as for the corresponding frames with nominal din>en 
sions. Henceforth, the frames with actual and nominal dimen- 
sions will be referred to as the actual and nominal frames, 
respectively. For the actual frames, spline interpolations 
through measured dimensions were used to generate the wall 
thicknesses and coordinates of the nodal points. Isoparametric 
finite elements were used to approximate the variations in 
stiffnesses and geometry. The one -dimensional models consid- 
ered herein are for the frames with nominal dimensions. The 
grids used for both the one-dimensional and two-dimensional 
models are described subsequently. 

T wo -dimensional models . An 18x8 grid was used for 
modeling the whole l-section frame. In this grid two elements 
were used to model each of the web, top and bottom flange 
sections. The part of the skin adjacent to the bottom flange 
section was treated as part of the flange. One element was used 
to model each of the two parts of the skin section extending 
beyond the bottom flange (see Fig. I ). The middle surfaces of 
the top flange and the web were taken to be their reference 
surfaces. The middle surface of the combined bottom flange 
and skin was taken to be the reference surface. 

An 18x7 grid was used for modeling the whole J-section 
frame. The distribution of the elements was similar to the 
l-section frame. Only one element was used to model the lop 
flange section (see Fig. I ). 

Totally clamped and partially clamped support condi- 
tions Mere considered. For totally clamped supports, all the six 
generalized displacements were restrained 
(u,'=u 2 '=w '=<(>! *=0 2 .=<|y =<)). The partially clamped 

conditions were obtained from the totally clamped case by 
successively removing the restraints on one. as well as on 
combinations, of the displacement and rotation components. 

( hie -dimensional models . A uniform grid of 24 elements 
was used in modeling each of the l-section and J-section 
frames The principal sectorial properties of the cross section 
were evaluated using the Fortran program listed in Coyette.™ 

Identification of Modes and Estimation of the Error in the 
One- Dimensional Model Predictions 

The two-dimensional models can be used to: a) identify 
the in plane, out-of-plane and coupled modes, and b) estimate 
the error in the predictions of the one-dimensional models. 
This is accomplished through decomposing the complementary 
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shnin energy. V\ Fqs. 9. assoc in tec l with each vibration mode, 
into three components, U,. U 2 ami U, (see Table I). Hie first 
two components. II, and U 2 are associated with the in-plane 
and out-of-plane stress resultants, respectively. 71ie third 
component, U v is associated with the stress resultants which 
are particular to two-dimensional plates and shells (not present 
in one-dimensional beam models). The in-plane and out-of- 
plane modes correspond to the modes for which IJ,/U C and 
IU/U* are close to I, respectively. The strongly coupled modes 
correspond to nearly equal values of U,/U c and U 2 /lJ t: . The 
ratio UyU 4 is indicative of the error in the one-dimensional 
model predictions. 

It is also useful to partition the total complementary 
strain energy, associated with each mtxle. into tlnee com- 
ponents. U,|. U w , IJhf representing the contributions of the top 
Mange, web, and bottom flange (including the skin). 

Comparison of Experimental and Finite Element Results 

The results of the experimental and numerical studies are 
summarized in Figs. 5 through 9 and Table 2 for the I-section 
I t ame, and in Figs. 10 through 14 and Table 3 for the J-section 
frame. For the finite element model three cases are considered, 
namely, totally clamped edges (with both translational and 
rotational restraints), partially clamped edges with not 

restrained, and partially clamped edges (with ti, - in the flanges 
and not restrained). 

The maximum and minimum values of the frequencies 
obtained by the two dimensional finite element model 
(corresponding to the totally damped and partially clamped 
edges) are shown in Figs. 5(a) and 10(a) along with the ex- 
perimental frequencies (See also Tables 2 and 3). Note that 
the experimental frequencies associated with mtxle 9 of the 
I section, and of the J-section. respectively, are close in fre 
quency. Mcxles for these frequencies have very close nodal 
locations. Also, the 12th iikhIc of the I section (see Table 2) 
was missed in the experimental survey which is indicative of 
die difficulty of determining the high frequency mcxles. The 
fact that only one of the multiple experimental frequencies with 
close ncxlal locations (nxxle 9) is predicted by the finite ele- 
ment model may be attributed to imperfections in lamination 
and material properties: and/or to geometric nonlinearilies 
which were not incorporated into the finite element mcxlel. In 
I igs. 5(b) and 10(b) bar charts are given for the frequencies 
obtained by two-dimensional nxxlels of the actual and nominal 
frames along with those of the one-dimensional model. 

In Figs. 6 and 1 1 bar charts are given showing the two 
decompositions of the complementary strain energies, associ- 
ated with the different vibration modes, described in the 
preceding subsection. The ordinates in Figs. 6(a) and 11(a) 
represent the ratios of U j /U° , U 2 /U c and Ut/U c , and the 
ordinates in Figs. 6(b) and 11(b) represent the ratios 
U,f/U c , 11*71)*’ . U|,f/lJ c for each of the modes. 

The nxxle shapes associated with the first five ex- 


perimental and analytical frequencies are shown in Figs. 7 and 
12. Two views are shown for the deformations associated with 
each mode: side view and lop view. Also shown are the nodal 
lines of the w' displacement on the top and bottom flanges. As 
can be seen from Figs. 7 and 12, the deformation patterns 
associated with higher modes are fairly complex. As men- 
tioned previously, the only experimental nodal lines monitored 
are those assix iated with gross in-plane, and gross out-of-plane 
motions. Generally, good agreement between the finite ele 
ment and experimental nodal lines is observed in these cases 
Other nodal lines, associated with localized deformations are 
shown only for the finite element solutions. 

The sensitivities of the vibration frequencies to the Fiber 
orientation angles of the top flange, web, and bottom flange 
and skin are depicted in Figs. 8 and 13. The ordinates in Figs. 
8 and 1 3 represent the sensitivity derivatives with respect to the 
indicated fiber angles. Each of the sensitivity derivatives is 
normalized by dividing it by the corresponding frequency of 
vibration. The sensitivities of the vibration frequencies to the 
material parameters F, , F,. G, , and G r| are shown in Figs 0 
and 14. The ordinates in Figs. 9 and 14 represent the sensitiv- 
ity derivatives with respect to the indicated elastic moduli. 
Each of the sensitivity derivatives is divided by the coi resjxxul 
ing frequency and multiplied by the corresponding elastic 
modulus. The effects of boundary conditions on the frequen- 
cies obtained by the two-dimensional finite element models are 
shown in Tables 2 and 3. 

An examination of the experimental and finite element 
results (Figs. 5 to 14 and Tables 2 and 3) reveals: 

1. Reasonably good correlation is observed between, 
numerical simulation and experiment for the l-section frame 
(see Fig. 5(b)). The ratios of the first five experimental fre- 
quencies to the corresponding finite element ones ranged 
between 0.92 and 1.02 (see Table 2). For the J-section frame 
the correlation is not as good (Fig. 10(b)). The corresponding 
ratios for the first five frequencies were 0.87 to 1.04 (see Table 
3). 

2. Most of the experimental frequencies for the l-section 
frame and the J-seciion frame are between those for the totally 
and partially clamped supports (with both ^ and u,' in the 

flanges not rest mined). This is particularly true for the higher, 
modes. For some of the n>odes the experimental frequencies 
are closer to the partially clamped support case (e g., modes 10, 

1 1 and 12. see Fig. 5(b)). For the I and J-section frames the 
finite element model predicted only one of the multiple ex- 
perimental frequencies with close nodal lines (mode 9). The 
other experimental frequencies were between those for the 
totally and partially clamped supports (with both <t> 2 and u, ' in 

the flanges not restrained, see Fig. 10(b)). 

3. The lowest five frequencies obtained by the one- 
dimensional model are reasonably close to those obtained by 
the corresponding two-dimensional model. This is particularly 
true for the J-beam where the errors in the predictions of the 
one -dimensional mcxlel were well below 10% (see Fies Sort 
and 10(a)). 
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4. Identification of the modes as in-plane or out-of-plane 
can best he accomplished by examining the energy com- 
ponents, U|/U c aiulUi/U 1 . associated with the in-plane and 
out-of-plane forces, respectively (Figs. 6(a) and 1 1(a)). Also, 
the error to he expected when using one -dimensional thin- 
walled beams can be estimated by computing the ratio of the 
energy associated with the forces neglected in thin- walled 
hesims to the total energy. U,/U* (see T igs. 6(a) and I 1(a)). 

5. The coupling between in-plane and out of plane 
deformations is more pronounced in the J -section than in the 
l-section frame. As an example, the first twenty modes for the 
I section frame had either U|/U l orllj/U 1 > 0.75. On the 
other hand, only modes I to 4, 6, X and 10 in the J -section 
frame had Uj/U c orlJ^AJ 1 > 0.75. For the higher modes 
neither the ratio U|/U l norl^/U 1 was close to I (see Figs. 
6(a) and I 1(a)). 

6. For the l-section frame, the contributions to the total 
energy of the top and bottom flanges far exceeded that of the 
web for any given mode. The ratio of the strain energy in the 
web to the total strain energy was less than 0.20 for the first ten 
modes and less than 0.2X for the succeeding ten modes (see 
Fig. 6(b)). For the J-section frame the strain energy in the web 
approached 0.4 of the total energy in some of the modes (see 
Fig. I 1(b)). 

7. For the l-section frame, the strain energy of the top 
flange is the dominant energy in the in-plane deformation 
nitxles and the strain energy of the bottom Mange (including the 
skin) dominates for the out-of-plane deformation modes (see 
Fig. 6(b)). 

X. The vibrational response of both the l-section and 
J-section frames is very sensitive to restraining the u, displace- 
ments of the flanges (and skin). It is somewhat sensitive to the 
rotational restraint on <|> 2 ' (see Tables 2 and 3). However, it is 

insensitive to restraining the displacement components and 
w\ and the rotation <)>, . 

9. The vibrational response of the I section and J-section 
is more sensitive to variations in the +45°, -45° fiber angles of 
the top flange than to variations in the 0° or 90° fiber angles. 
The variations in the 0° and 90° fillers of the web and the 
bottom flange have a noticeable effect on some of the nnxles. 
but their effect is generally less than that of the 45 n , -45° l ibers 
(sec Figs. X ami 13). The vibrational response is also more 
sensitive to variations in the elastic nnxluli F, and (i, , than to 
any of the other material coefficients (see Figs. 9 and 14). 

10. I he sensitivity of the vibration frequencies with 
respect to variations in both F, and G, , is almost the same for 
all the modes (see Figs. 9 and 14). lliis may be attributed to 
the quasi-isotropic lamination used for both the flanges and the 
web. It suggests the feasibility of replacing the quasi -isotropic 
composite, in the one-dimensional thin-walled beam model, by 
an equivalent isotropic material, as was done in the present 
study. 


Comments on Sources of Errors 
Sources of Errors 

The detenni nation of natural frequencies and modes 
from vibration tests and numerical models involves numerous 
possible sources of discrepancies or errors which are related to 
mechanical and equipment limitations as well as to theoretical 
and physical assumptions. The errors in vibration tests include 
inexact equipment calibration, excessive noise, manufacturing 
variations, incorrect transducer locations and operation in a 
region of nonlinearity of the response. Numerical nxxleling 
errors can be attributed to inaccuracies in estimated material 
properties and insufficient nxxleling detail. In the present 
study care was exercised in collecting and recording (he 
vibration test data, and in the selection of the numerical model. 
However, nominal material properties and layups (fiber orien- 
tation of the different layers) were used in the numerical model. 
The sensitivity analysis helped in identifying the material and 
lamination parameters that need to be accurately determined. 

Concluding Remarks 

A detailed study was made of the effects of variations in 
lamination and material parameters on their vibrational charac- 
teristics of thin-walled composite frames. Ilie structures 
considered are semicircular, thin- walled frames with I and J 
cross sections. The Manges, web and skin of the stiffeners have 
quasi-isotropic laminations with filler orientation being combi- 
nations of ± 45°, 0° and 90° layers. Two computational nxxlels 
are used for predicting the vibrational characteristics. In the 
first model, the Manges and webs of the stiffeners were 
modeled by using two-dimensional shell (and plate) finite 
elements. The second model was a finite element discretization 
of the one-dimensional Vlasov's-type thin-walled beam theory. 
A mixed formulation was used with the fundamental unknowns 
consisting of both the generalized displacements and stress 
resultants (or internal forces) in the frame. The frequencies and 
mtxles predicted by the computational models are compared 
with those obtained from experiments. A detailed study was 
made of the sensitivity of the vibrational response to variations 
in the fiber orientation, material properties of the individual 
layers, and boundary conditions. On the basis of this study the 
following conclusions are justified: 

1. For some of the higher vibration modes the ex- 
perimental frequencies for thin-walled frames are generally 
between those for the totally and partially clamped supports. 

2. Identification of the modes as in plane or out-of -plane 
can best be accomplished by examining the energy components 
associated with the in-plane and out-of-plane forces. Also, the 
minimum error to be expected when using one-dimensional 
thin-walled beams can be estimated by computing the ratio of 
the energy associated with the forces neglected in thin-walled 
beams to the total energy. 

3 For quasi -isotropic composite frames the vibration 
frequencies, associated with the lower modes, can be accurately 
predicted by isotropic one dimensional beam model (with 
effective elastic moduli). The accuracy of predictions is 


5 



dependent on the cross- -sectional distortions during the beam 
deformations. As the cross-sectional distortions increase, the 
degradation of accuracy becomes more pronounced. 

4. The vibrational response of thin- walled semicircular 
frames is very sensitive to restraining the displacement compo- 
nent of the flanges along the length of the frame. It is some- 
what sensitive to the restraint on the associated rotational 
component. However, it is less sensitive to restraining the 
other displacement and rotation components. 

-5- The vibrational response of thin-walled composite 
Irames with quasi isotropic laminations is more sensitive to 
variations in the +45", -45" liber angles of the top flange than to 
variations in the 0‘ or 90" fiber angles. Variations in the O' and 
90° filers of the web and the bottom flange have a noticeable 
effect on some of the modes, but their effect is generally less 
than that of the 45 . -45 fibers. The vibrational response is 
also more sensitive to variations in the material coefficients E, 
and ( ij | than to all other coefficients. 

b. 1 he sensitivity of the vibration frequencies with 
respect to variations in both E, and G, , is almost the same for 

all the nxxles. This may be attributed to the quasi- isotropic 
lamination used for both the flanges and the web. It suggests 
the feasibility of replacing the quasi-isotropic composite by an 
equivalent isotropic material in the one-dimensional thin 
walled beam analysis, as was done in the present study. 
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one-dimensional beam nuxlel 
bending stress resultants in the two 
dimensional nuxlel 

consistent and generalized mass matrices 
for an individual element (see Eqs. 3) 
extensional stress resultants in the two- 
dimensional model 


N, 

. <i, 
R 

R| 


t r 




Hi 


U. V, w 


axial force in the one dimensional beam 
nude! 

transverse shear forces in the one- 
dimensional beam model 
transverse shear stress resultants in the 
two-dimensional model 
radius of curvature of the centerline of the 
frame (used in one dimensional beam 
model)' 

outer radius of curvature of the frame (see 
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displacement components of the two 
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coordinate directions 

displacement components of the two- 
dimensional model in the x, ,x 2 ,x 3 ' 

coordinate directions 

vector of nodal displacements 
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Table I.- Decomposition of total complementary energy, U c 
into components. 
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Global coordinate 
system 



E^ = 2.0 x 107 psi 
E t = 1.7 x 10® psi 
G^j = 9.3 x 10® psi 
Gjj = 6.51 x 10® psi 


J LT 


= 0.38 


P = 0.058 lb/in.3 
Nominal layer thickness = 0.005 in. 

Fiber Orientation 


NL = 8 


[+45/0/90] 


NL = 16 : [±45/0/90] 2s 


Figure 1 Thin-walled composite frames and coordinate systems used in present study. 
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Web 


Figure 2.- Sign convention for generalized displacements and stress resultants in two-dimensional model. 



Figure 3.- Sign convention for generalized displacements and stress resultants for one-dimensional model. 
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Figure 4.- Photograph of thin-walled semi-circular 

graphite-epoxy specimens and experimental 
equipment. 



(a) Two-dimensional and one-dimensional 
beam model results. 


Frequency, 

Hz 



(b) Experimental and bounding two-dimensional 
model results. 

Figure 5.- Comparison of finite element and experimental 
frequencies for the thin-walled composite 
frame with I cross-section. 
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(b) U tf , U w , and components. 

Figure 6. - Energy components in the different vibration 
modes of the thin-walled composite frame 
with I cross-section. 
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Top View 


Side View 



c<>2 = 31 .9 (29.7) Hz 


Nodal Lines on Top Experimental Node Lines 

and Bottom Flanges (Bottom Flange) 



o» 3 = 37.5 (35.9) Hz 





Figure 7.- Mode shapes associated with the lowest five frequencies for the thin-walled 

composite frame with I cross-section. Parentheses are experimental frequencies. 
Others are clamped finite element model results. 
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(a) Top flange. 



(b) Web. 



(c) Bottom flange and skin. 


Figure 8.- Sensitivity of vibration frequencies to fiber 

orientation angle in the flanges and web of the 
thin-walled composite frame with I cross- 
section. 



Figure 9.- Sensitivity of vibration frequencies to variations 
in material characteristics of the thin-walled 
composite frame with I cross-section. 
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Modto 

(a) Two-dimensional and one-dimensional 
beam model results. 

Frequency, Hz 



(b) Experimental and bounding two-dimensional 
model results. 


Figure 10.- Comparison of finite element and 

experimental frequencies for the thin-walled 
composite frame with J cross-section. 
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(a) U 1 , U 2 , and U 3 comoponents. (b) U tfl U w , and U bf components. 


Figure 1 1 Energy components in the different vibration modes of the thin-walled 
composite frame with J cross-section. 


Top View Side View Nodal Lines on Top Experimental Node Lines 

and Rnttnm Flannac ( Bottom Flanno\ 



(o 2 = 36.9 (32.1) Hz 



0 ) 4 = 79.2 (69.0) Hz 



« 5 = 91 .4 (79.0) Hz 


Figure 12 - Mode shapes associated with the five lowest frequencies of the thin-walled 
composite frame with J cross-section. Parentheses are experimental 
frequencies. Others are clamped finite element model results. 
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(a) Top flange. 
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Figure 14.- Sensitivity of vibration frequencies to 

variations in the material characteristics of the 
thin-walled composite frame with J cross- 
section. 



(c) Bottom flange and skin. 

Figure 13.- Sensitivity of vibration frequencies to 

variations in the fiber orientation angle in the 
flanges and web of the thin-walled 
composite frame with J cross-section. 
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